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1 Introduction 

It may seem surprising that a Hamiltonian lattice theory whose single-site states only 
have the quantum numbers of mesons and baryons can be equivalent to a theory of 
free massless quarks, but it is true. I will show that this result follows directly from 
the application of the Contractor Renormalization Group (CORE) method[|l| to an 
appropriately chosen free fermion theory. 

The original motivation for the computation I will present was the desire to apply 
CORE to lattice quantum chromodynamics (QCD) and show that the picture which 
emerged from older strong-coupling calculations [|| also applies to the weak coupling 
regime. The new feature of CORE which allows this question to be dealt with non- 
perturbatively is that CORE, as opposed to earlier Hamiltonian real-space renormal- 
ization group approaches, allows one to retain only gauge-invariant block-states (i.e., 
states in which no flux leaves a block) and still define a non-trivial renormalization 
group transformation. 

Fluxless states are of particular importance to the strong coupling limit of a lat- 
tice gauge-theory because every link which carries nonvanishing flux contributes an 
energy proportional to g^, so that fluxless states have the lowest energy. Moreover, if 
a lattice theory allows for single-site color-singlet states (e.g., theories which exhibit 
explicit fermion doubling, Wilson fermions, domain-wall fermions, or theories based 
upon variants of the SLAC-type of derivative) the number of degenerate fluxless states 
is huge and in the limit g^ —>■ oo these states will all have zero energy. In the case 
of lattice QCD, single-site meson states (i.e., color-singlet quark anti-quark states), 
single-site baryon states, and single-site multi-meson and baryon states consistent 
with the exclusion principle, are all of this type. The huge degeneracy among these 
fluxless states is lifted in order l/g"^ and, for a nearest-neighbor derivative, perturb- 
ing in l/g"^ leads to an effective Hamiltonian which has the form of a generalized 
Heisenberg antiferromagnet. An immediate consequence of this result is that chiral 
SU{3) X SU{3) is spontaneously broken in the groundstate of this theory; another gen- 
eral result is that the theory has an approximate SU (12) symmetry which is broken 
if one adds next-to-nearest neighbor terms to the fermion derivative. 

Although it is attractive to rewrite strong-coupling QCD in terms of states which 
have the quantum numbers of mesons and baryons, establishing the relevance of the 
chiral-symmetry prediction and the approximate SU{12) results to the continuum, 
small g'^, limit is problematic. CORE allows one to systematically study this ques- 
tion by truncating the Hilbert space to this set of strong-coupling states, obtaining a 
renormalized Hamiltonian which is valid for all values of g. The important question 
which must be answered is whether truncation to this set of states biases the com- 
putation and incorrectly forces the strong coupling results of confinement and chiral 
symmetry breaking. One way to show that this is not the case is to apply the same 
truncation scheme to free fermion theory and show that it leads to a renormalized 
Hamiltonian which has the same physics as the global color-singlet sector of the free 
theory. While the general theorem on which CORE is based guarantees this result 



will hold if the retained states have a non-trivial overlap with the relevant low-lying 
states of the free-field theory, a-priori nothing forces this overlap to be nonvanishing; 
establishing this fact requires a calculation. 

This paper does the requisite calculation for a theory of free massless fermions in 
1 + 1 space-time dimensions with a nearest neighbor fermion derivative. It will be 
clear that the extension of this calculation to higher dimensions and other derivatives 
is straightforward. While the nearest-neighbor theory exhibits species doubling and 
is anomaly free, a fact which makes it useless for studying the physics of theories such 
as the Schwinger model, it is simple to use and is relevant to the question of whether 
or not truncating to the natural strong-coupling states makes it impossible to obtain 
the correct weak-coupling physics. The explicit calculation shows that things work 
as expected for the nearest-neighbor theory and that the renormalized Hamiltonian 
takes the form of a generalized frustrated antiferromagnet which, perforce, has all 
the physics of the charge zero sector of the original free fermion theory. 

Since there is no substantive difference between the physics of free relativistic 
lattice fermions and the physics of this unusual frustrated antiferromagnetic system 
(which doesn't, at the short distance level, have any relativistic fermions) it is inter- 
esting to ask which theory is fundamental ? Clearly, at this level there is no way to 
decide the issue. It will be apparent from the calculation that the same mapping can 
be carried out for different choices of fermion derivative with similar results. In other 
words, the couplings will vary in strength but the general form of the renormalized 
Hamiltonian will be the same. From a renormalization group point of view this says 
that there is a surface in coupling constant space of the generalized frustrated anti- 
ferromagnetic system where all Hamiltonians flow to the same free massless fermion 
fixed point. Furthermore, it will be clear that turning on gauge fields in the origi- 
nal problem produces the same sort of Hamiltonian with different coefficients. Thus, 
the gauge theory is also hidden inside this same system or, in other words, all of 
these different theories are different phases of the same generic Hamiltonian. As a 
class, frustrated antiferromagnets are systems which have recently come to be of some 
interest because of their possible connection to high- Tc superconductors, 0, so this 
identification of a more general class of HAF's (for a specific couplings) to a theory 
of free relativistic fermions, or relativistic fermions interacting through a gauge-field, 
has interest beyond its application to QCD. 

2 A Brief Review of CORE 

The CORE method consists of two parts, a theorem which defines the Hamilto- 
nian analog of Wilson's exact renormalization group transformation and a set of ap- 
proximation procedures which render nonperturbative calculation of the renormalized 
Hamiltonian doable. 

CORE replaces the Lagrangian notion of integrating out degrees of freedom by 
that of throwing away Hilbert space states. This is accomplished by defining a pro- 



jection operator, P, which acts on the original Hilbert space, 7i and whose image is 
defined to be the space of retained states Tirct = P'H. The foundation of the CORE 
approach is a formula which relates the original Hamiltonian, H, to the renormalized 
Hamiltonian which has, in a sense which was made precise in Ref. [|l|, exactly the 
same low energy physics as H. This equation is 

W'^^ = hm [[T(t)2]]-V2 [[T{t)HT{t)]] [[T(t)^]]-l/^ (1) 

t — >OQ 

where T(t) = e~*^ and where [[O]] = POP for any operator O which acts on Ti. 
A similar formula can be written to define the renormalized version of any other 
extensive operator. The new, renormalized, operator is guaranteed to have the same 
matrix elements between eigenstates of H'^'^'^ that the original operator had between 
eigenstates of H. 

2.1 The cluster expansion 

Generally one cannot evaluate Eq.|I] exactly, however it is possible to nonperturba- 
tively approximate the infinite lattice version of H'^'^'^ to any desired degree of accu- 
racy. This works because W'^'^, as defined in Eq.|l|, is an extensive operator and has 
the general form 

oo 
^ren^^^conn^^-^^^) (2) 

r=l 

where each term, h™"^^{j,r), stands for a set of range-r connected operators based at 
site j, all of which can be evaluated to high accuracy using finite size lattices. The 
explicit definition of the connected range-r operator, /i™°'^(j, r), depends upon the 
details of the truncation procedure. In what follows I will limit discussion to the case 
of a one-dimensional lattice, since this is what I will need to discuss the free-fermion 
theory. A detailed discussion of the general methodology can be found in Ref. [|I| . 

2.2 The Approximation Procedure 

Three steps define the nonperturbative approximation scheme for computing H'^'^^: 
first the truncation procedure; second the subtraction procedure used to convert 
the evaluation of Eq.|l] on finite sublattices to the operators /i™°'^(p, r); third the 
method for evaluating the t —>■ oo limit in Eq.|l|, without explicitly computing either 
[[T^(t)]]~^/^ or [[T{t)HT{t)]]. I will heuristically review each of these steps in turn. 

First some notation. In what follows I deal with one-dimensional spatial lattices 
whose sites are labelled hj —oo < j < oo. 1 assume that there are A^-states corre- 
sponding to each site j of the lattice which I denote by \(pa{j))j, where a = 1 . . . A^. 
A basis for the full Hilbert space H will be generated by taking tensor products of 
these A^-states per site over all sites j. 

CORE allows a wide choice of truncation procedures, however I will limit myself to 
one which appears to work well in a large number of cases. First, divide the lattice into 



disjoint blocks Bp each having ub sites and keep a small number states per block. The 
way to choose which states to keep is to diagonalize the block-Hamiltonian (i.e., that 
Hamiltonian obtained by restricting H to only those terms which are contained within 
any one of the blocks Bp) and throw away all but its M lowest lying eigenstates, where 
M < A^"^. If we let Hbp denote the block Hamiltonian and |\E'q,(p)) for a = 1..M its 
eigenstates, then the projection operator P is 



P = IlPp 

p 

M<N-^B 

Pp = Y. I*«(P))(*«(P)I (3) 

Given P it only remains to compute W^^. 

Generally the lattice on which the renormalized Hamiltonian is defined is thinner 
than the original lattice in that each site p on the new lattice corresponds to a block of 
sites Bp on the old lattice, however this need not be the case. In the free-fermion case, 
to be discussed in the next section, we will thin the states associated with a single 
site of the original lattice and map the original theory into an equivalent theory with 
the same number of sites but with a Hamiltonian which has a very different form. 

To define the cluster expansion begin by defining the range- 1 term in H 
/i™°'^(p, 1), to be 



ren 



h'"^^{p,l) = PHb,P 



M 



E^«I*"(p))(*"(p)IKp (4) 



\a=l 



P P 

^ p -^ p 



where HBp\'^a{p)) = -Sa|^«(p)) and where 

P = 
Pp = X{Pi (5) 

The range-2 connected operator /i'^°°'^(p, 2) is defined by subtracting /i™°'^(p, 1) and 
^conn^p _l_ -^^ ]^^ from H^^^Bp, -Bp+i), the operator obtained by evaluating the t -^ oo 
of Eq.|l] for H restricted to the two adjacent blocks {Bp, -Bp+i}. Note, in this case the 
notation [[O]] stands for P = PpPp+iOPpPp-^i. The explicit definition of h'^"'^^{p, 2) is 

/i=°-(p, 2) = P^^^^, H^'\Bp, Bp+,) Pp-^p+, - /.--(p, 1) - /.--(p +1,1). (6) 

where, in analogy to Eq.^ I define 

Pp,p+i= n Pi (7) 



Similarly, the range three operator /i'^°°'^(p, 3) would be obtained from the following 
formula 

-/i'=°""(p + 2, 1) - /i^°"°(p, 2) - /i^™^°(p + 1, 2) (8) 

In this case one must subtract the three different ways of embedding the connected 
range- 1 computation and the two-different ways of embedding the connected range-2 
computation in the three block problem. (As before, H^^\Bp, -Bp+i, -Bp+2) is the oper- 
ator which results when one restricts H to the three adjacent blocks {Bp, -Bp+i, -Bp+2} 
and then evaluates Eq.|I[) 

Although it is possible to numerically evaluate Eq.|l] for any multi-block sublattice 
and extract the limit t ^ 00 by taking large vales of t (see for example Ref.|]I|) this 
is not necessary The fact is that each term in the cluster expansion can be computed 
from a knowledge of the tensor product states (which span the space of retained states 
for the multi-block problem and the eigenvalues and eigenstates of the corresponding 
multi-block Hamiltonian. A general proof of this assertion appears in Ref . |]l| , however 
the basic ideas are summarized in the following theorem. 

Theorem: Let H^ be a single block Hamiltonian and let P be the projection 
operator which corresponds to keeping its M lowest lying eigenstates |\E'a); 
furthermore, let H denote the Hamiltonian of an r-block sublattice and let 
the M^ tensor product states formed from the states |\E'q,) span the space of 
retained states. Then the t — > cxd limit of the equation which defines the 
renormalized multi-block Hamiltonian can always be written as 

H'''' = lim[[T(t)2]]-i/2r[T(t)i7T(t)]][[T(t)2]]-^/2 

t — ^■oo 

= RH^i^^R^ (9) 

where R is an M^ x M^ orthogonal matrix, i?"!" its inverse and H^iag is a diagonal 
matrix whose entries are the eigenvalues of those M^ lowest lying eigenstates 
of H which appear in the expansion of the retained states. 

To clarify what is meant by the iVf lowest lying eigenstates which appear in the 
expansion of the retained states 1 will consider two simple examples. The first example 
corresponds to the simplest truncation procedure one can imagine, i.e. choosing 
M = 1 and truncating to a single state. In this case the theorem is trivial since ikf = 1 
and so the space of retained states for that multi-block system is one dimensional. 
The fact that R and -R^ are orthogonal matrices means their single matrix element 
must be 1 and so, as long as the single retained state has an overlap with the ground- 
state, -ffdiag must simply be the groundstate energy of the multi-block Hamiltonian. 
To prove this assertion it suffices to rewrite Eq.|I] as 



mt)) 



(^|e-2t^|^) 



H"'"'' = lim(^(t)|iJ|^(t)) (10) 



and then expand |\l/(t)) in a complete set of eigenstates of H. 

The second example, M = 2 and r = 2, exhibits all essential features of the 
general case. Obviously, in this case iVf = 4 and so [[T^{t)]]-^/^ and [[T{t)HT{t)]] 
are 4x4 matrices, each of which becomes singular in the limit t -^ oo, although 
their product is well defined. To understand why the product is well defined and has 
the form shown in Eq.^ it is convenient to expand the four retained states in terms 
of exact eigenstates of the two-block problem and on the basis of this expansion, 
construct an orthogonal transformation R which renders the evaluation of the limit 
in Eq.p straightforward. If we denote the four retained states as I^E'i), 1^2), I ^3) and 
|\E'4) we can write their expansion in terms of block eigenstates as 

1^1) = ao|0o) + ai|0i) + 02102) + ••• 

1^2) = &o|0o)+&l|0l)+&2|02) + --- 
1*3) = Co|0o)+Cl|0l)+C2|02) + ... 
1^4) = rfo|0o)+rfl|0l)+C?2|02) + -.- (11) 

where the states |0„) correspond to eigenstates of the block Hamiltonian with energies 
e„. Assume that the states are arranged in the order of increasing energy, so that 0o 
is the groundstate of the block Hamiltonian, |0i) the first excited state, etc. 
Applying T(t) to each of these states we obtain 



T(t) ^1) 


= aoe-*^« |0o) + a,e~''' 0i) + 026"*^^ ^s) + • • • 


T(t) ^2) 


= 6oe"*'° 0o) + e-*^^6i |0i) + fese"*^^ |02) + • ■ ■ 


T(t) M/3) 


= Coe-*^° 0o) + cie-*^^ |0i) + Cse"*^^ 02) + • • • 


T{t) ^4) 


= rfoe-*^« |0o) + die-''' 0i) + ^26"*^^ |02) + • • • 



(12) 

The reason it is convenient to make an orthogonal transformation on the states |\l'j) 
is that in the t —* 00 limit those states in Eq.[l2| for which the coefficient of |0o) is 
nonvanishing will, up to a normalization factor, contract onto the same state I^q), 
rendering [[T(t) i7T(t)]] and [[T^it)]] singular. By multiplying [[T{t)HT{t)]] by the 
factors of [[T(t)^]]~^'^ we correct for this situation, but it is not at all obvious why or 
how this works in the original tensor product basis. 

To avoid this problem with the large t limit we change basis, defining states |xi) to 
1X4), which are orthonormal linear combinations of the states |\E'i) to |^E'4), having the 
property that each state contracts onto a distinct eigenstate of the block Hamiltonian; 
i.e., 

IXl) = ao|0o) +«l|0l) +«2|02) +"3103) + ••• 

1X2) = A|0l)+/?2|02)+/53|03) + ... 

1X3) = 72I02) +73I03) + ••■ 

1X4) = 53|03) + .-- (13) 



Applying T(t) to these basis states yields 



Wlxi) = 


aoe ' 


^^«|0o) + aie-*^^|0i) + a^e^'^' 02) + age-'^^l^s) + 


{t)\X2) = 




Ae-*^^ 0i) + /52e-*^^|02) + /53e-*^^|03) + • 


Wlxs) = 




72e-*^^02)+73e-*^^03) + . 


WIX4) = 




Sse-'-' 03) + . 



(14) 

from which it is clear that, up to a normalization factor, each of the states \xa) 
contracts onto a different eigenstate of the multi-block Hamiltonian. Furthermore, 
we have 

H T{t) \xi) = «oeoe-*^" |0o) + ftieie"*^! |0i) + «2e2e-*^^ |02) + a3e3e-*^^ |03) + • • • 
HT{t) 1x2) = /3ieie-*^i|0i) + /32e2e-*^l02) + Psese-'^' Ih) + ■■■ 

HT{t)\x3) = 72626-*^^ 102) +73636-*^^ 103) + ... 

HT{t)\x,) = 53e3e-*^l03) + . .(15) 

Given Eq.|l^ and Eq.^ one can easily analyze the t -^ oo behavior of Eq.^ To get 
a feeling for the way in which this works let us simplify the discussion by, for the 
moment, assuming that [[T^(t)]] and [[T (t) HT (t)]] are 2x2 matrices obtained by 
sandwiching these operators between the states |xi) and \x2)- In this case we have 



[[T'im 



mmrm 









In general [[T^(t)]] is a matrix of scalar products and its eigenvalues are guaranteed 
to be positive, so it is guaranteed to have an inverse square root. On can explicitly 
construct the inverse square root by writing [[T^(t)]] as 

[[T\m = U{t)^D{t)U{t) (16) 

where D{t) is a diagonal matrix whose entries are the eigenvalues of [[T'^(i)]] and U{t) 
is the matrix whose columns are the normalized eigenvectors corresponding to those 
eigenvalues. Given this decomposition 

[[T2(t)]]-l/2 ^ ^(^)t /^(t)-l/2 ^(i) (17) 

Note, since D(t) is diagonal, D{t)^^/'^ is also a diagonal matrix whose entries are the 
inverse square roots of the corresponding entries in D{t). 

Fortunately, since all we really need is the behavior of the product in Eq.^ as t gets 
large, we don't have to do all this work. It suffices to define a t-dependent rescaling 

8 



of Ixi) and \x2) which guarantees that each state converges, as t ^ oo, to the lowest 
lying eigenstate of the block Hamiltonian which appears in its expansion in terms of 
multi-block eigenstates. More specifically, multiplying |xi) by the factor e'^''*/<^o and 
1X2) by ei t/ I3i yields the result 

Xi) = l0o) + -e-(^-^»)*|0i) + ... 



'At)) - 


= nt) — \xi) 


At)) - 


peit 

= nt) — \x2) 
pi 



I01) + ... 

from which it follows that in the limit of large t, 



[[T\t)]] 



mt)HTit)]] 



1 + ... 4e"('i"'o)* 



-'o 



^ g-(<:i-eo)t I I 

ao 



eo + • • • 



£1 g g-(<:i-<:o) t 




which establishes the theorem for this 2x2 case. It should be clear that the same 
sort of rescaling of \xi) to Ixa) will establish the result for the real 4x4 case. It 
is important to reiterate that the construction the matrix of eigenvalues D and the 
construction of the orthogonal transformation R is done directly from a knowledge 
of the eigenvalues of the block Hamiltonian and the expansion of the tensor product 
states in eigenstates of the block Hamiltonian; at no point is it necessary to deal 
with [[T(t)]] for finite values of t. This means that when dealing with large blocks 
and many eigenstates, techniques such as the Lanczos method, which finds only the 
relevant lowest lying eigenstates starting from the tensor product states, can greatly 
reduce the computer resources needed to solve the problem. 

3 Free-Field Preliminaries 

To set up the computation presented in the next section requires some notation. 
The system under discussion is a 1 + 1-dimensional Hamiltonian lattice theory. The 
fermion field is taken to be a two component operator ipa{j) with a = 1,2. The 
Hamiltonian has the form 

H=j:l^'in-J2)^i^3^n (19) 

where 5 is a general hopping term having the property 

S'{ji-J2) = -S'{j2-Ji). (20) 



^U) = ^.(J) + 4(' 



and o"3 is the 2x2 matrix 

I -I) (-) 

As is customary I introduce creation and annihilation operators bj,bj,dj,dj and 
define the nearest neighbor derivative, 6 (j), as follows: 

5'{1) = -5\-l) = l 

6'ij) = ifjVl (22) 

Making these substitutions the Hamiltonian takes the form 

^ = E ^(^Vi^. - b%^i) - lid]+idj - d]dj+i) (23) 

In terms of these operators the total electric and axial-charge are defined by 

Q = Y.i%-d]dj) 
j 

Q' = T.ibh + d]d,-l) (24) 

j 

For each site j, there are two electrically neutral states. The first is the state |0j) 
which is annihilated by bj and dj. The other chargeless (i.e. locally gauge- invariant) 
state is \±j) = bjdj\Oj). The projection operator onto the space of retained states is 
defined in terms of these states to by the product 



(25) 





p = HP. 




P, = 0,)(0, + ±,)(±, . 


Fourier transforming 




bk 


i i 


dk 




the Hamiltonian in Eq.|23| becomes 


H 


= 1. r dksm{k){bibk-didk) 

ZTT J-tt 



(26) 



(27) 

It follows from Eq.^ that the ground-state of the theory is the state obtained by 
filling the negative energy sea; i.e., 

|vac) = n&U4|0) (28) 

fc>0 
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where |0) is the state annihilated by all the b^s and dkS ( or equivalently, the 6/s 
and dj's). 

The only additional free-fermion formula which I need has to do with diagonalizing 
Eq.^for a finite lattice where j = 1..N. As with all such quadratic Haniiltonians one 
only has to diagonalize the N x N matrix Kj^j^ = i/2 6 (ji — J2), which can be done 
exactly for any value of A^. If A^ is even the eigenvalues of K look like the eigenvalues 
of the infinite volume problem in that 

e{k) = sm{kp) (29) 

where kp = 7rp/2{N + 1) for p any odd integer such that — A^ < P < A^- For 
the case of odd A^, the eigenvalues are once given by Eq.^ however in this case 
the values of kp are given by kp = 7Tp/{N + 1) where now p is any integer such 
that -(A^ + l)/2 < p < {N + l)/2. Note that in both cases the range of k is 
— 7r/2 < k < tt/2 in distinction to the case of a periodic or infinite lattice. Although 
the issue of how doubling works in the finite volume open boundary condition case is 
interesting, time does not permit going into it here. Suffice it to say that even though 
the eigenvalues are given by Eq.^ and undoubled, nevertheless the theory has both 
left and right movers for each component of tp and is not a chiral theory. The last 
piece of information needed in order to carry out the full computation is the formula 
for the eigenvectors of K, which for all values of A^ are given by 

UkM = / / , ie-^'^' - i-iye'""') (30) 

V2(ivT2) 

A final fact concerning the nearest-neighbor free-fermion theory is that it possesses 
a global SU{2) symmetry whose generators are given by 



ba 


= > ] 2^1^3^"^3V^: 


s+ 


j 


S- 




S3 


= hh% + d]d,- 



J 



(31) 

This symmetry is very useful for checking the results of the necessary finite volume 
calculations. 



4 The Computation 

To compute the connected range-1 term iJ'^°°'^(j, 1) one diagonalizes the single-site 
Hamiltonian and selects a set of retained states from among its lowest lying eigen- 
states. Since there are no range-1 terms in Eq.^ the single-site Hamiltonian is zero 

11 



and so, independent of the choice of retained states /i™'^°(j, 1) = 0. As noted earher, 
I will retain the locally gauge-invariant states 

|0,); |±,) = 6]4|0,). (32) 

which form a spin-1/2 doublet with respect to the global SU(2) defined in Eq.pT|. 

The first non-trivial connected operator which contributes to H'^'^'^ is the range-2 
operator /i'^°°'^(j, 2). The four retained states associated with a connected two-site 
block are 

|00); |±0); |0,±); |±±) (33) 

which decompose into a spin-1 and spin-0 multiplet under the global SU{2). Since 
each of the states of definite spin and third-component of spin must contract onto 
unique eigenstates of the two-site problem, changing from the tensor product states 
to the spin basis amounts to finding the rotation matrix R. The combinations of 
definite total spin and 3-component of spin are: 



sL = i 


^3 = 


-1 


|00) 


sL = 1 


^3 = 





^(l±o)- o±)) 


sL = i 


^3 = 


1 


±±) 


SL = 1 


^3 = 





^(l±o) + |o±)) 



(34) 

Note the unusual minus sign in the spin-1, 3^ = state, comes from the alternating 
minus sign which appears in the definition of S+ and S'_ in Eqs.^ and Eq.pT|. 

To complete the computation of /i'^°'^"(j, 2) we only need to find the eigenvalues of 
the full two-site Hamiltonian which correspond to the lowest lying spin-0 and spin- 
1 states. Since the Hamiltonian, Eq.|23|, only has terms which absorb a particle or 
antiparticle at site j and create one at j ± 1 it follows that |00), which is annihilated 
by the absorption operators, has zero energy. The same is true of the state |±±), 
which is not annihilated by the absorption operators but which has no room to move a 
particle or antipartile to a neighboring site. It follows from the SU{2) symmetry that 
the state -^(|±0) — |0±)) also has zero energy. The only eigenvalue to be determined 
is that associated with the lowest spin-0 state in the sector which has one particle and 
one antiparticle. For the case of two sites the two allowed values for the momenta 
are k = ±7r/6 and e{k) = ±sin(7r/6) = ±1/2; thus the lowest energy charge zero 
eigenstate of the two-site problem is 6_^/g(iLg|0, 0), having an energy E2 = —1. It 
follows from this that, in the total spin basis 



Hrii,i + l)=R 



f 0\ 





V -I J 
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R^ (35) 



where R is the rotation matrix which transformed the original tensor products into 
the total spin basis. This 4x4 matrix can always be rewritten in terms of the sixteen 
matrices Ma/3 = cra{i)(J/3{i + 1) where a, /5 = ... 3 and ctq is the unit matrix and 0"i, 
(72 and (Js are the Pauli spin matrices. In particular if we write 

iJr(^,^ + l)= E Ca/3M,^ (36) 

«,/3=0 

then it follows from the S'f/(2)-symmetry and Eq.|3^ that, in the original tensor prod- 
uct basis, H^'^ii^i + 1) has the form 

/i^°'^°(2, i) = Hr{i, z + 1) = -ll + s{i) ■ s{i + 1) (37) 

since h'^°'^^{l,i) = 0. Thus, up to range-2 terms, the full renormalized Hamiltonian is 

H^^^^ = -^V + Y.m-s{^ + ^) (38) 

i 

Actually, this form of the renormalized Hamiltonian gives the correct spectrum, but it 
is not what one obtains directly by taking traces of the form tT{sx{i)sx{i + ^)H2^'^{i, i + 
1)) as indicated in Eq.|36|. The nonvanishing traces are 

^tT{a,{^)a,{^ + l)Hr{^,^ + l)) = -\ 

h.T{a■,{^)a■S + m^{h^ + ^)) = \ (39) 

and so the form of if|'^'^ directly obtained from our definitions has the form 



—^ + E(-^i(0 ^i(^ + 1) - ^2(0 S2[i + 1) + s^{i) s^{% + 1) 

i 



(40) 



Note, however, that this Hamiltonian can be brought into the form of a Heisenberg 
antiferromagnet by the rotation 

= ne~*"''^''^'^ (41) 

which, for every odd lattice site maps Si(2j + 1) ^ — si(2j + 1) and S2(2j + 1) — > 
— S2(2j + 1) but leaves S3(2j + 1) unchanged. In the calculations which follow we will 
see that, given the phases I have chosen for for the states |0j) and |±j), the terms in 
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H^'^^ will generically have the form shown in Eq.^ and will be brought into standard 
form by application of O. 

While the computation of /i™°'^(l, i) and h'^°'^^{2, i) are very simple, the calculation 
of h'^°^'^{n,i) for n > 3 requires more work. Fortunately, it is easy to automate this 
task using MapleV and the computation of all terms out to and including h'^°'^^{5,i) 
takes less than two minutes on a desktop computer. I will now describe the process 
for the case /i™°'^(3,'i), since the general computation proceeds along the same lines. 

The allowed momenta for a three-site sublattice are k = — 7r/4, 0, 7r/4; the corre- 
sponding particle energies are sin(— 7r/4) = — 1/v^, sin(O) = 0,sin(7r/4) = 1/v^; and 
the antiparticle energies are — sin(— 7r/4) = 1/v^, — sin(O) = 0, — sin(7r/4) = — 1/a/2. 
As the Hamiltonian contains only hopping terms, the number of particles and an- 
tiparticles are separately conserved; thus, of the eight possible tensor product states 
(which define the set of retained states) the states |000) and |ib ± ±) are eigenstates of 
the three-site Hamiltonian of energy zero. Furthermore, of the six-remaining states, 
the three two-particle states |±00), |0 ± 0) and |00±) and the three four-particle 
states |± ± 0), |±0±) and |0 ± ±) can be treated separately, since when expanded in 
terms of a complete set of eigenstates for the three-site problem they have no states 
in common. 

The nine possible particle-antiparticle eigenstates of the three site problem are 

\k,l) = bl4\000) (42) 

with eigenenergy e{k,l) = sm{k) — sin(/). Expanding the three retained states in 
terms of these eigenstates is straightforward since 

6J4|000) =Y.Uk{j)ui{j)bld]\000) (43) 

k,l 

thus, the coefficient aki of the state bid} is just 

aj{kJ)=J2Mj)ui{j) (44) 

k,l 

Given these formulae we can now discuss how to compute the matrices i?,ifdiag 
and R^ from the overlap matrix, which gives the expansion of the retained states in 
terms of the eigenstates of the three-site Hamiltonian. The eigenenergies of the one 
particle states ordered by energy are: 

e(-7r/4,7r/4) = -v^ , e(-7r/4, 0) = -1/v^ , e(0, 7r/4) = -1/v^ 
e(-7r/4, -7r/4) = , e(0,0) = , e(7r/4, 7r/4) = (45) 

e(7r/4,0) = l/v^ , e(0, -7r/4) = -1/^2 , e{n/A,-7r/A) = ^2 

The first overlap matrix giving the expansion of the retained states in terms of these 
eigenstates is 

O^ = (0 ± 0| I .50 -.5 -.5 I (46) 



(±0 


/ .25 


.35 


.35 


.25 


.25 


.5 


.35 


(0±0| 


.50 








-.5 


-.5 








(0 0± 


V.25 


-.35 


-.35 


.25 


.25 


.5 


-.35 
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Focusing on the first column of this matrix it is easy to construct a 3 x 3-rotation ma- 
trix, Ri which takes the first column into a vector with one nonvanishing component; 

i.e. 

/ .61 -.20 -.20 .41 \ 

RiO^=\ -.39 -.39 -.37 -.37 -.38 -.39 (47) 

V -.32 -.32 .45 .45 .45 -.32 / 

The fact that the first column has one nonvanishing entry means that it is the only 
state which has an overlap with the three-site groundstate, which is the first objective 
we wished to achieve. Now, the fact that the last two entries in the second column are 
nonvanishing means that we are still in a situation where both the second and third 
rotated retained states contract onto the first excited state of the three-site problem. 
Clearly we can now perform a rotation R2 which is the identity on the first rotated 
state and rotates the second two states so that in the new rotated basis only the 
second element of the second column oi R2R2OV is nonvanishing; i.e.. 



61 








-.20 


-.20 


.41 








.5 


.5 











.5 











-.58 


-.58 


-.58 






R^R^O^ = .5 .5 .5 (48) 

V -.58 -.58 -.58 / 

Now we have finally achieved our original goal, each of the rotated states will contract 
onto a unique eigenstate of the three-site Hamiltonian; i.e., the first rotated state 
contracts onto the state h_^ ud\ u\Q{i{i) , the second state contracts onto {b^^udl + 

6^)4/4) |000)/v^ and the third state contracts onto (6^^^/4^1^/4+6^4+^1/44/4) |000)/v^- 
Thus, in this new rotated basis the diagonal Hamiltonian in the three-site one pair 
sector is 

/ -^2 \ 
i7D(l-pair)= -I/V2 (49) 

V y 

Given the explicit form of the rotation matrix R = R2R1 we can compute the Hamil- 
tonian in the original tensor product basis to be 



i73^*'^(l-pair) = RHD(l-pair)Rt = -.471 -.943 -.471 (50) 

V .118 -.471 -.589 J 

A similar calculation can be done for the two-pair sector of the three-site problem. 
Here the three retained states are |± ± 0), |±0±) and |0 ± ±). Obviously, as in the 
one-pair case, these states are of the generic form 

\ji,j2) = bldlbldl\000) (51) 

and can be directly rewritten in terms of the eigenstates of the three-site Hamiltonian 
using the operators 6^ and d] as 

\ji,J2) = J2 UkMui,{ji)uk2{J2)ui2{J2)bi^dlbi^dl\000) 

ki,k2,li,l2 
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ki<k2 li<l2 

Z J2i'^hiJl)Uh{J2)-Ul,{ji)ui,{J2))\ki,k2jlj2) (52) 

ki<k2 h<h 

where by |fci, A;2, ^i, ^2) we denote the state bk-^di_^bl^dijOOO) for the particular ordering 
fci < ^2 and /i < /2- With this choice of a complete set of basis vectors it is easy 
to read off the overlap matrix; it is not surprising, given the SU{2) symmetry of the 
nearest neighbor problem, that it is identical to the overlap matrix of the one-pair 
sector and so the same procedure can be applied to bring it to the desired form. 

Note that the no-pair state |000) and the three-pair state |± ± ±) are each exact 
eigenstates of the three-site problem with eigenvalue zero so 

Hl''''{0 - pair)(i, i + 1, i + 2) = ff3^"^(3 - pair)(i, i + 1, i + 2) = (53) 

Given these results it is now trivial to expand Hl^'^ as 

Hrii, t + l,t + 2) = J2 c^p^i^M^A^ + ^>S + 2) (54) 

0/37 

by taking traces; thus, for example, the coefficient of the identity operator is 

cooo = ^tr(a„(z)(T^(z + \)a^[% + 2)iJr(z, z + 1, ^ + 2)). (55) 

o 

Proceeding in this way we obtain, for the case of three sites. 

To convert this to a connected contribution we have to subtract the two ways of 
embedding the two-site problem into the three-site problem, obtaining 



d— sitesconn J— sites cy 

'^000 ~ '^000 ^"^00 



,3— sitesconn ^3— sites o 2— sitesconn 

00 ^"^00 

.53033 ... - 2(-l/4) = -.03033 . . . (57) 



which shows that the range-3 connected contribution to the ground-state energy den- 
sity is quite small in comparison to the range-2 contribution. Similarly we can com- 
pute the coefficients of the other operators which can appear in Eq.^ Thus, we 
have 

3— sites 3— sites 3— sites 9Qc;7n 

•-110 ~ —'-220 — '^330 — .ZOd/U... 

3— sites 3— sites 3— sites OQc;7n 

•-Oil — —"-022 — '-033 — .zoom... 

4-f- = =4-f^^ = 4-f^^ = -.05893... (58) 

and all other possible Cq,/3^ vanish. 

As we already noted, for nearest neighbor spin-spin interactions, which is what 
Clio, C220 and C330 parameterize, there should be a difference in sign between C330 and 
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both Clio and C220, which is the case; the same is true for cqh, C022 and C033. Note 
however, the next to nearest neighbor coefficients Cioi, C202 and C303 should have the 
same sign and they do. 

To obtain the connected coefficients we recall that 

Hr^'i^, t+i,t+2) = i7r(^, ^ + 1, ^ + 2) - //r°(^, ^ + 1) - i/2°"°(^ +1,^+2). (59) 



Using the explicit form for H^'^^i, i + 1, i + 2), H', 
and collecting like terms we obtain 



conn 

2 y-i 



i,i + l) and/J2'°°"( 



i + lA + 2) 



110 

conn 

220 

conn 

'330 

conn 

Oil 

conn 

'022 

conn 

'033 

conn 

101 

conn 

202 

conn 

303 



3— sites 
110 

3— sites 
220 
3— sites 
330 

,3— sites 
Oil 
3— sites 
'-022 
3— sites 



C 



C 



C 



c; 



c, 



033 



2— sites 
' '^11 

2— sites 
" '-22 

2— sites 
' '^33 

2— sites 
' '^11 

2— sites 
" '-22 

2— sites 
" '-33 



0.014297.. 

0.014297.. 

-0.014297, 

0.014297.. 

0.014297.. 

-0.014297, 



0.05892... 
0.05892... 
0.05892... 



(60) 



all other coefficients vanish. 

Thus, recalling the for a = 1, 2, 3 s^ = o'a/2, adding up the range-2 and range-3 
connected coefficients and rotating Si and S2 to minus themselves on alternate sites, 
we can rewrite the range-3 renormalized Hamiltonian as 

^ren ^ _o.28033\/ + ^ [0.9428s(i) ■ s{i + 1) + .23570s(i) ■ s{i + 2)] (61) 



which has the advertised form of a frustrated Heisenberg antiferro magnet. The impor- 
tant things to notice about the result of this calculation is that the range-3 corrections 
to the coefficients of operators which appeared at range-2 are small and that the co- 
efficients of the new operators which appear for the ffist time at range-3 are typically 
smaller than the ones which appear at range-2. 

This process can be carried out in the same way to compute W^™^. The result 
out to and including range-5 contributions is 



H' 



_.31099y + ^iJ2-body(0 



+ E^' 



4-body(^j 



H- 



2-body(?j 



H. 



4-body 



= 0.80001 s(i 

= 0.03559 s(i 

- 0.08033 s(i 

+ 0.03403 s(i 

+ 0.02595 s(i 



s{i + 1) + 0.23492 s{i) ■ s{i + 2) - 0.01915 s{i) ■ s{i + 3) 
s{i + 1) s{i + 2) ■ s{i + 3) 
s{i + 2) s{i + 1) • s{i + 3) 
s{i + 3) s{i + 1) ■ s{i + 2) 
s{i + 1) s{i + 2) ■ s{i + 4) 
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+ 0.00339 s(i) 

- 0.01159 s(i) 
+ 0.05189 s'(i) 

- 0.03289 s(i) 

- 0.01159 s(i) 

- 0.00732 s(i) 

- 0.03251 s(i) 

- 0.00732 sfi) 



s{i + 1) s{i + 3 
s{i + 2) s{i + 1 
s{i + 2) s{i + 3 
s{i + 3) s{i + 1 
s{i + 3) s{i + 2 
s{i + 4) s{i + 1 
s{i + 4) s{i + 1 
s{i + 4) s{i + 2 



s(i + 4) 
s{i + 4) 
s{i + 4) 
s(i + 4) 
s{i + 4) 
s(i + 2) 
s(i + 3) 
sii + 3) 



(62) 



The important fact to notice about all of these calculations is that at each stage 
the tensor product states contract onto the lowest energy states of the correspond- 
ing block Hamiltonian, thus proving that, if one computes out to infinite range, the 
renormalized Hamiltonian will describe zero-charge sector of the free system. Fur- 
thermore, note that if one uses the same method to calculate the two-point fermion 
antifermion correlation function CORE will guarantee that its renormalized operator 
has the same matrix elements in the frustrated HAF ground-state that the original 
operators had in the free-fermion vacuum state. 



5 Discussion 

Although, in principle, there are an infinite number of terms in H'^^^ Eq.|6^ shows that 
the coefficients of the longer range two-body terms, as well as the coefficients of the 
four-body and higher terms drop off quickly; thus, as advertised in the introduction, 
we see that the zero-charge sector of the free (doubled) fermion theory maps into 
a generalized frustrated antiferromagnet. Frustrated, in this context, simply means 
that the coefficient of the range-two and range-three one-body terms are both positive 
and the range-four term is negative. To understand what is going on physically con- 
sider a Neel state in which the spins are oriented |. . . |||| . . .). It is clear, restricting 
attention to the nearest-neighbor s{i) ■ s{i + 1) terms, that the fact that the coeffi- 
cient of these terms is positive means that this orientation of the spins minimize this 
contribution to the energy. This, of course, implies that spins separated by two sites 
should point in the same direction. However, this maximizes instead of minimizes 
the range-three spin-spin term since its coefficient is also positive. Note also the the 
range-four term also works against the range-two term since it wants spins separated 
by three sites to point in the same direction, which is opposite to what one obtains 
from miminizing the range-two term. The presence of the four-body (and higher 
body) operators further confuses the issue. Since little is known about the physics of 
frustrated antiferromagnets, it is of some interest to know that this generalized class 
of theories are related to free fermion theories. In fact, I will now argue that this is 
true for a non-trivial range of parameters. 



To begin note that it is easy to modify the free fermion derivative by adding 
hopping terms which jump three sites, five sites, seven sites, etc., in such a way that 
the infinite volume kinetic energy has the general form 

R 

e{k) = ^ flr sin((2r + l)k) (63) 

m=l 

This form of e{k) is still symmetric about k = it/2 and exhibits the same sort of 
fermion doubling as is exhibited by the nearest neighbor term. In this case it is a 
little more difficult to carry out the computations of the previous section, but they 
go through in essentially the same way, except that the coefficients of the operators 
in iJ'''='i are changed. It therefore follows that there is more than one generalized 
frustrated antiferromagnet which is equivalent to a theory of free fermions. Since all 
of the various fermion derivatives can be defined so that the slope of e{k) is unity 
for small values of k, it follows that all of these different theories must flow, under 
real-space renormalization group transformations, to a single flxed-point generalized 
frustrated antiferromagnet. This means that there is a nontrivial surface in the pa- 
rameter space of these generalized frustrated antiferromagnets which at low energies 
describe the same system of free massless fermions. 

Having said this, it is important to point out that not all of these frustrated 
antiferromagnets correspond to theories of free fermions. As I pointed out in the 
introduction the same mapping can be carried out for the case of an interacting gauge 
theory. In fact, the original motivation of this work was to do exactly this for lattice 
QCD and analyze the weak coupling regime. Obviously, constructing the mapping 
which takes the interacting gauge theory into this class of generalized frustrated 
antiferromagnets is technically more challenging, but it must exist. From this it 
follows that for some choice of couplings the generalized frustrated antiferromagnet 
corresponds not to free fermions, but rather to a theory with fermions and gauge 
flelds, neither of which is apparent when one looks at the theory at the single-site or 
few-site level. 

A final issue which needs to be addressed is how things work for higher dimension 
and/or non-abelian gauge theories. Clearly, for the case of doubled free fermion theo- 
ries in 3 + 1-dimensions there are no new technical challenges; however, the number of 
single-site gauge-invariant states which can be constructed is larger since the particles 
and antiparticles come with two different spins. As in the 1 + 1-dimensional case, 
these states form an irreducible representation of a symmetry group, in this case a 
six-dimensional representation of SU{4) and truncating to this set of states, using 
CORE to compute the renormalized Hamiltonian, leads to a generalized frustrated 
SU{4) antiferromagnet (i.e., one where the s ■ s terms are replaced hj Q ■ Q terms, 
where the Q's are now the representation matrices for the six-dimensional representa- 
tion of the generators of SU{4)). In this case using generalized to modify the phrase 
frustrated antiferromagnet is even more appropriate since in this theory a new kind of 
interaction term appears in the renormalized Hamiltonian. This happens because the 
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six-dimensional representation of SU{A) is not the fundamental representation and 
so operators of the form {Q ■ QY can appear, and these operators can significantly 
alter the landscape of possible phases which the theory can have. A discussion of how 
such terms can modify the behavior of an antiferromagnetic theory will be given in a 
forthcoming paper on the Haldane conjecture.^ Note, generally all terms permitted 
by the global symmetries of the problem can and will appear. 

As to the case of free quarks in 3 + 1-dimensions, things are again different in 
detail, but the general results are similar. Here we are interested in all the color 
singlet states which can be formed from quarks and anti-quarks on a single site. 
If one restricts to just mesons, i.e. all states formed from color-singlet quark anti- 
quark states, then the states on a single site fall into an irreducible representation of 
SU {12), which is a symmetry of the doubled theory. Once again, if one truncates to 
just this set of states the resulting renormalized Hamiltonian must be a generalized 
frustrated SU{12) antiferromagnet. Once again generalized antiferromagnet must be 
taken in the broadest possible sense since new terms of the form discussed for the 
case of a single fermion in 3 + 1-dimensions will appear, but now one will also get 
new manybody terms related to the possibility of having Casimir operators beyond 
the simplest quadratic Casimir operator. Thus one should expect that theories which 
correspond to quarks in interaction with color-gauge fields to have a very rich struc- 
ture. Obviously, all remarks concerning the option of choosing different kinetic terms 
for the lattice fermion theory apply equally to the free massless quark case, and so 
there must be a whole surface in the space of couplings for the generalized frustrated 
theory which all flow to the same fixed point. 

As for real lattice QCD with quarks interacting with gluons, it has been shown 
that at strong coupling the single-site colorless states are all degenerate to leading 
order in g"^ and that this enormous degeneracy is lifted in order 1/g'^, see Ref.0]. 
Working to this order in the doubled theory one obtains, exactly as in the case of 
the l+!-dimensional theory, a generalized antiferromagnet. Thus, we see that for 
this truncation algorithm the free field theory and the strong coupling theory can be 
obtained from one another by varying parameters in the general renormalized Hamil- 
tonian. This tells us that the full generalized frustrated SU {12) antiferromagnet 
theory has regions in which the degrees of freedom are no longer free massless quarks, 
but are instead quarks interacting with color gauge fields. Note, adding the baryonic 
states which can be made from qqq configurations leads to a further generalization of 
the antiferromagnetic system in which there is more than one irreducible represen- 
tation of the symmetry group for each site in the lattice and correspondingly more 
complicated interactions. All of these couplings become significant as one moves to 
weak coupling but in the limit of very large g^ the states with baryons in them split 
away from the states with only mesons in order 1/ g"^. 

Finally, returning to the point raised in the introduction; namely, although this 
exercise was originally motivated by the desire to show that CORE provides frame- 
work within which the old strong coupling treatment of QCD can be extended to the 
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weak coupling regime, at this point the question arises as to whether or not one can 
tell if in fact quarks and gauge-fields are really fundamental. More thought has to be 
given to the question of whether or not one can meaningfully distinguish lattice QCD 
from the frustrated antiferromagnet. Clearly much work remains to be done to settle 
how various features of the various undoubled theories, such as the anomaly, work in 
practice. Also, it remains to be seen how all of these undoubled theories differ from 
one another. 
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